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Vanadium dioxide (VO2) with its unique sharp resistivity change at the metal-insulator transition (MIT) has 
been extensively considered for the near-future terahertz/infrared devices and energy harvesting systems. 
Controlling the epitaxial quality and microstructures of vanadium dioxide thin films and understanding the 
metal-insulator transition behaviors are therefore critical to novel device development. The metal-insulator 
transition behaviors of the epitaxial vanadium dioxide thin films deposited on AI2O3 (0001) substrates were 
systematically studied by characterizing the temperature dependency of both Raman spectrum and Fourier 
transform infrared spectroscopy. Our findings on the correlation between the nucleation dynamics of 
intermediate monoclinic (M2) phase with microstructures will open a new avenue for the design and 
integration of advanced heterostructures with controllable multifunctionalities for sensing and imaging 
system applications. 



Vanadium dioxide (VO2) undergoes a metal-insulator transition (MIT) with the resistivity changing sharply 
in 4-5 orders of magnitude just above the room temperature (341 K). Such unique characteristic can be 
advantageously utilized in microelectronic devices, terahertz devices, and energy harvesting systems, 
etc^"^. The MIT is considered to be accompanied by a structural phase transition from monoclinic to rutile. 
Recently, the formation of a metastable monoclinic phase of VO2 during the phase transition has attracted much 
attention, not only from the point of view of the theoretical understanding of the MIT mechanism, but also from 
the technological consideration in control of the transition properties for device applications. The most stable 
monoclinic phase of VO2 at room temperature is referred to as the Ml phase, with the j5-angle of the unit cell as 
122.6°. In this structure, the V atoms are paired and twisted across the octahedral edgesl The localized electron 
state in this structure is considered to induce the insulating state. However, during the MIT of the insulating Ml 
phase to the metallic rutile phase, another metastable monoclinic phase with the j5-angle of 91.88°, named as the 
M2 phase, has been reported. This metastable M2 phase is normally induced by either doping or strain. In the M2 
structure, half of the V atoms pair but do not twist as in the Ml phase. In contrast, the other half of the vanadium 
atoms twist, but remain equidistant along the resulting zigzag chain^"^. 

To understand the nature of the M2 phase, both theoretical and experimental efforts have been initiated. It was 
proposed that the uniaxial stress along [110]r direction would stabilize the M2 phase^. Experimentally, investi- 
gations on the M1-M2 transition were mainly focused on single crystalline one- dimensional nanostructures 
because it is easier to apply a mechanical stress on these structures^'^"^^. There are only a few reports on the 
observation of the M2 phase in the VO2 thin films, although stress is widely considered as an important feature in 
thin film systems. Okimura et at. conducted in situ X-ray diffraction (XRD) studies on the epitaxial VO2 films on 
AI2O3 (0001) grown by pulsed laser deposition, indicating the formation of the M2 phase during the heating 
process^^. More recently, they also found that the M2 phase could be stabilized in polycrystalline thin films at 
room temperature by using a sputtering technique with inductively coupled plasma^^'^^ These results have clearly 
demonstrated that thin film VO2 would be a good model system to study the intermediate phase. However, due to 
the very complex microstructures of VO2 thin films and related phase transitions, fundamental understanding of 
the correlations between the microstructures and the occurrence of the M2 phase in VO2 thin films is still lacking. 
In this paper, we report our effort to systematically study the occurrence of the M2 phase in epitaxial VO2 thin 
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Figure 1 | (a) XRD 6 — 26 scans for samples Th96 and Thl56, respectively; 
(b) XRD scans for samples Th90 and Thl56, respectively. 

films, using both Raman and Fourier transform infrared (FTIR) 
spectroscopies at various temperatures. To fully understand the 
mechanisms of the formation of the M2 phase, the high resolution 
X-ray diffraction (HRXRD) technique was employed to analyze the 
strain states and the microstructures. 

Results 

Briefly, the VO2 thin films were grown on AI2O3 (0001) substrates by 
using a polymer- assisted deposition (PAD) technique as reported 
earlier^^"^^. Different thicknesses of VO2 films were synthesized: 
one set of the sample (labeled as sample Th90) had a thickness of 
90 nm and the other set of the sample (labeled as sample Thl50) had 
a thickness of 150 nm. Figs. 1 (a) and (b) show the XRD patterns 
from 0 — 20 and (\) scans on the films, respectively. The results 
indicate that both samples are crystallized in monoclinic VO2 phase 
(Ml) at room temperature and epitaxially grown on the substrates. 
The epitaxial relationship between the thin film and the substrate was 
determined to be VO2 (001)||Al2O3 (0001) and VO2 [010]||Al2O3 
[11 20]. 

Raman spectroscopy was carried out at various temperatures to 
check the occurrence of the M2 phase. Samples were heated from 300 
to 370 K during the measurement. Each data point was the average 
from 10 repeated records with the integral time of 10 s. The variation 
of the temperature was controlled to be less than 0.5 K during the 
measurement. According to the reported results the fingerprint 
for the M1-M2 phase transition is the blueshift of the Raman peak 
from 618 cm"^ (the Ml phase) to 642 cm"^ (the M2 phase). The 
blueshift is resulted from the change of Ag vibrate mode in VO2 unit 
cell. As shown in Fig 2(a), both samples (Th90 and Thl50) exhibit 
almost the same spectra at room temperature. The spectra are 
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Figure 2 | (a) Raman spectra for samples Th90 and Thl50 at room 
temperature, respectively; (b) Raman spectra for sample Th90 at 
temperatures ranging from 300 to 360 K; (c) Raman spectra for sample 
Thl50 at temperatures ranging from 300 to 360 K. 

identical to those from the pure Ml phase in the bulk single crystal- 
line V02^^'^°. This result is reasonable considering the good crystal- 
line quality of the as-prepared epitaxial films. With the increase of the 
temperature, as shown in Figs. 2(b) and (c), the peaks in the Raman 
spectrum become weaker and finally disappear. This corresponds to 
the first order transition of crystalline structure from the low sym- 
metry of monoclinic to the high symmetry of rutile. If we look at the 
detailed Raman spectra during the temperature increase, it is inter- 
esting to notice the difference between Th90 and Thl50 at around 
334 K. Specifically, occurrence of the M2 phase of VO2 with the 
Raman peak shifting from 618 cm"^ to 642 cm"^ could be clearly 
observed in sample Thl50 instead of sample Th90, as shown in 
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Figs. 2(b) and (c). To verify the results, the Raman spectra between 
520 and 720 cm"^ for both samples measured at various tempera- 
tures were carefully fit with Lorentz functions. Figs. 3(a) and (b) show 
the fitting results for Th90 and Thl50, respectively. The peaks cen- 
tered at 577, 618, 642 cm"^ correspond to the characteristic Raman 
peaks from the Eg mode of sapphire, and the Ag modes of the Ml and 
M2 phases of VO2, respectively. It is found that both samples are 
actually composed of mixing phases of Ml and M2 during the MIT. 
However, the M2 peak started to occur at a much lower temperature 
for sample Thl50 (305 K) than Th90 (332 K). Moreover, the M2 
peak is comparable or even higher than the Ml peak for Thl50 when 
the temperature is higher than 330 K. On the contrary, for Th90, the 
Ml peak is dominant throughout the whole measuring range of 



temperatures. In other words, both the amount of the M2 phase 
and the temperature window for the occurrence of the M2 phase 
are larger in Thl50 than Th90 during the MIT. 

To have a better understanding of the phase transition, IR trans- 
mission properties of the samples were measured by using the IR 
light with different wavelengths at different temperatures. Both the 
as-prepared samples, Th90 and Thl50, exhibit strong modulation of 
IR transmission at room temperature in the wavelength range from 
2500 to 5000 nm, as shown in Fig. 4(a). The transmittance becomes 
lower as the film thickness increases, possibly due to the stronger 
absorption from a thicker film. Fig. 4(b) is the IR transmittance 
measured at various temperatures with a wavelength of 4000 nm. 
Hysteresis loops were observed from both samples as the results of 
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Figure 3 | Lorentz fitting results for the Raman spectra in the range from 520 cm"^ to 720 cm~^ measured at various temperatures for samples Th90 (a) 
and Thl50 (b), respectively. (Black dots are the original data, black lines are the Lorentz fitting line, purple lines indicate the Eg mode of sapphire, red lines 
indicate the Ag mode of the Ml phase of VO2, and blue lines indicate the Ag mode of the M2 phase of VO2.) 
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Figure 4 | (a) FTIR spectra for samples Th90 and Th 150 at 313 K and 
353 K, respectively; (b) IR transmittance versus temperature for samples 
Th90 and Thl50, respectively; (c) differential curves of transmittance 
versus temperature for samples Th90 and Thl50, respectively (the dotted 
lines indicated the fitting results of these differential curves using Gaussian 
functions). 

the intrinsic MIT properties of VO2. The differential curves of the IR 
transmittance versus temperature are plotted in Fig. 4(c). If we fit 
these differential curves using Gaussian functions, as shown in the 
dotted lines in Fig. 4(c), two peaks located at different temperatures 
(332.1 K and 334.5 K) are observed in sample Thl50 during the 
cooling process, whereas no such obvious split is observed in sample 
Th90. The two peaks in sample Thl50 are most likely to be related to 
the M1-M2 transition and M2-rutile transition, respectively, which 
is in good agreement with the M2 occurrence observed in the Raman 
spectra. 

To understand the mechanism of the occurrence of M2 phases, the 
strain states/distributions of the thin films were characterized by 



using HRXRD. Similar to the strain distribution analysis for the other 
epitaxial oxide thin films^^"^^, the strain states can be determined by 
characterizing the changes in related diffraction peak positions, 
which reflect the changes in the lattice parameters. The reciprocal 
space maps (RSMs) were recorded around AI2O3 (0006) and VO2 
(002) for samples Th90 and Thl50, as shown in Figs. 5(a) and (b), 
respectively. It is found that the positions of VO2 (002) for both 
samples are almost the same, i.e., at about 3400 arcseconds lower 
than the diffraction spots from AI2O3 (0006). The lattice parameters 
of the c-axis for both samples (Th90 and Thl50) are estimated to be 
5.3730 A. Gomparing to the lattice parameter of bulk VO2, 5.3820 A 
(JGPDS 82-0661), a compressive strain of —0.17% was found along 
the c-axis for both samples. Under the assumption of constant unit- 
cell volume, a tensile strain along [011]mi is expected. By tilting the 
samples with the setting of % = 44.8416°, XRD 0 — 20 scanning 
patterns were collected to measure the in-plane lattice parameters 
along the [011]mi direction of the as-prepared VO2 thin films, as 
shown in Fig. 5(c). The [011]mi lattice parameters were determined 
to be 3.208 A and 3.209 A for samples Th90 and Thl50, respectively. 
Gomparing with its bulk parameter (3.2034 A), tensile strains of 
0.14% and 0.17% for samples Th90 and Thl50 were estimated along 
the [011]mi direction, respectively. 

Discussion 

The [01 1]mi would be changed to [1 10] r when the monoclinic phase 
Ml is changed to rutile structure^^. Therefore, the tensile strains may 
induce the formation of the M2 phase in both samples, Th90 and 
Thl50^\ Nevertheless, the temperature window for the occurrence of 
the M2 phase can be quite narrow due to the very weak strains. 
Although the values of the tensile strain in samples Th90 (0.14%) 
and Thl50 (0.17%) are very close to each other, the occurrence 
window for the M2 phase is very different. More exactly, the amount 
of M2 phase and the temperature window for the occurrence of the 
M2 phase are much larger in Thl50 than Th90, leading to the clear 
direct observation of M2 phase in sample Thl50 rather than in Th90 
under the exactly same experimental set up, as demonstrated by the 
temperature dependency of Raman spectrum and FTIR spectrum 
shown above. 

To clarify the different occurrence window of the M2 phase in 
these two epitaxial VO2 thin films, peak spreading of the as-grown 
films were analyzed by HRXRD co and <\> scans. As shown in Figs. 5(a) 
and (b), the full width at half maximum (FWHM) of 00 — 20 scans 
around VO2 (002) (along they-axis in Figs. 5(a) and (b)) are almost 
the same for samples Th90 and Thl50, indicating the similar global 
strain states of both films. However, the FWHM of the 00 scan (along 
the X-axis in Figs. 5(a) and (b)) for the Thl50 sample is much wider 
than that for the Th90 sample. Similar results (as shown in Fig. 5(d)) 
were observed in the 00 scan rocking curve of VO2 (Oil) where the 
sample was tilted by x = 44.8416°. Theoretically, the FWHM of the (D 
scan (Pm) comes from several sources^^: 



(1) 



where Po is the intrinsic width of X-ray, Pd is from optical path of the 
X-ray equipment, Pj. depends on the radius of curvature of the sam- 
ple. Pi depends on the grain size of the sample, and P^ is the broa- 
dened width due to the screw dislocations of the sample. Under the 
same measurement set-up, equation (1) clearly suggests that the 
difference in the FWHM of the 00 scan should be from the broadened 
width induced by screw dislocations or grain boundaries, i.e. the 
value of (Pi + Pa) in sample Thl50 is higher than that in sample 
Th90. On the other hand, it is known that the density of edge dis- 
locations can be estimated from the peak broadening of the (j) scan. 
As shown in Fig. 5(e), the larger FWHM of the (j) scan peaks for 
sample Thl50 than that for sample Th90 suggests higher density of 
edge dislocations in the sample Thl50. The higher densities of both 
screw dislocations and edge dislocations may induce higher density 
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Figure 5 | (a) RSM around AI2O3 (0006) and VO2 (002) for sample Th90; (b) RSM around AI2O3 (0006) and VO2 (002) for sample Thl50; (c) tilted XRD 
6 — 26 scans {% = 44.8416°) for samples Th90 and Thl50, respectively; (d) Gaussian fittings of co rocking curves of VO2 (Oil) for samples Th90 and 
Thl50, respectively; (e) Gaussian fittings of <^ scans of VO2 (Oil) for samples Th90 and Thl50, respectively. 



of domain/grain boundaries in sample Thl50. These domain/grain 
boundaries may result in a strong impact on the nucleation of the M2 
phase^^'^°'^^'^^. The low energy at the domain/grain boundaries would 
promote the nucleation of the new phase. Thus, the richer domain/ 
grain boundaries in sample Thl50 would increase the opportunities 
for the formation of the M2 phase, which either increases the amount 
of M2 domains or broadens the temperature windows for the occur- 
rence of an M2 phase. This explains the clear direct observation of the 
M2 phase during the MIT in sample Thl50 rather than in sample 
Th90. 

In conclusion, the epitaxial VO2 thin films were prepared on 
AI2O3 (0001) using a polymer- assisted deposition technique. The 
intermediate M2 phase was clearly observed in the thicker thin films 
by temperature dependent Raman and FTIR spectroscopies. 
Although HRXRD demonstrated that both the as -prepared samples 
show very similar tensile strains along the [110]r direction, the 
higher density of dislocations and domain/grain boundaries may 
account for the easier nucleation of the M2 phase in the thicker film 
(150 nm thick). These results are helpful to understand and control 
the occurrence of an M2 phase as well as the resistivity and IR switch- 
ing properties during MIT, which are critical to the device applica- 
tions of VO2 thin films. 

Methods 

The VO2 thin films were grown on AI2O3 (0001) substrates by using a polymer- 
assisted deposition (PAD) technique. 2.0 g polyethyleneimine (PEI, average Mn ~ 
60000, Mw « 750000), 3.5 g ethylene diamine tetraacetic acid (EDTA) and 1.4 g 
NH4VO3 were dissolved in 40 mL distilled water to form a homogeneous polymer 
precursor solution. The solution was then purified in an Amicon filtration (Amicon 
8050) unit and concentrated to get a precursor solution containing 12.53 mg V 
element per mL, as measured by an Inductively Coupled Plasma Optical Emission 
Spectrometer (ICP-OES). The as-prepared solution was spin-coated on the AI2O3 
(0001) substrates. After that, the precursor films were thermally treated at 480°C in 
the reduced atmosphere of forming gas (H2 2% and N2 98%) for 1 hour. The 



spin-coating and thermal treatment process was repeated for various times to achieve 
films with different thicknesses. 

The XRD characterization was conducted using a Bede Dl diffractometer with Cu 
Kqj radiation (/I = 1.54060 A). HRXRD measurements were set by adopting a Si (220) 
dual channel collimator and monochromator system to provide a highly collimated 
Cu Kqji beam with a wavelength spread less than the natural line width of Cu K^i and a 
divergence of 12 arcseconds. The detector acceptance aperture was strictly limited by 
a precision slit with the width of 50 |j,m. 

Raman spectra were recorded at various temperatures using a Renishaw Micro- 
Raman Spectroscopy System. The wavelength of the Raman exciting laser used in this 
measurement was 532 nm. 

The infrared ray (IR) transmittance versus temperature was measured by using a 
Fourier transform infrared spectrometer (FTIR-8400S) equipped with a home-made 
temperature controlling unit. 
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